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[i]  Discovered  almost  four  decades  ago,  the  midnight 
temperature  maximum  (MTM)  with  typical  magnitudes  of 
50-100  K  has  been  regularly  observed  by  satellite  and 
ground-based  instruments  in  the  tropical  upper 
thermosphere.  Although  several  mechanisms  have  been 
suggested  to  explain  the  phenomenon,  previous  attempts  to 
reproduce  it  with  comprehensive  thennosphere-ionosphere 
models  have  been  unsuccessful.  First  long-term  simulations 
with  the  Whole  Atmosphere  Model  (WAM)  reveal  the 
presence  of  a  realistically  prominent  MTM  and  reproduce 
the  salient  features  of  its  daily,  seasonal,  and  latitudinal 
variability.  Preliminary  analysis  indicates  that  the  feature 
may  be  traced  down  to  the  lower  thermosphere,  where  it  is 
manifested  primarily  in  the  form  of  an  upward  propagating 
terdiumal  tidal  wave.  Its  spectrum  expands  to  higher-order 
zonal  wavenumbers  and  frequencies  and  its  phase  advances 
to  near  midnight  higher  up  in  the  thermosphere.  Possible 
mechanisms  generating  this  wave  may  involve  nonlinear 
interactions  between  other  tidal  harmonics  originating  in  the 
middle  and  lower  atmosphere.  Our  results  thus  suggest  that 
the  MTM  is  yet  another  phenomenon  driven  by  dynamical 
links  between  the  lower  and  upper  atmosphere  and 
ionosphere.  Citation:  Akmaev,  R.  A.,  F.  Wu,  T.  J. 
Fuller-Rowell,  and  H.  Wang  (2009),  Midnight  temperature 
maximum  (MTM)  in  Whole  Atmosphere  Model  (WAM) 
simulations,  Geophys.  Res.  Lett.,  36,  L07108,  doi:10.1029/ 
2009GL037759. 

1.  Introduction 

[2]  The  MTM  was  first  reported  in  radar  measurements 
of  ion  temperature  in  the  thermosphere  over  Arecibo  (18°N) 
in  the  early  1970s  [Harper,  1973],  In  1977-1978  system¬ 
atic  observations  were  made  at  low  latitudes  equatorward  of 
±20°  by  the  Neutral  Atmosphere  Temperature  Experiment 
(NATE)  instrument  on  the  Atmosphere  Explorer-E  (AE-E) 
satellite  [Spencer  et  al.,  1979],  An  AE-E  climatology 
describing  seasonal  and  horizontal  variations  of  the  MTM 
was  compiled  in  an  altitude  range  of  250-370  km  for  this 
period  of  low  solar  activity  [Herrero  and  Spencer,  1982],  It 
is  generally  consistent  with  ground-based  optical  observa¬ 
tions,  primarily  in  the  tropical  American  sector,  often 
indicating  a  poleward  propagation  of  the  feature  with 
typical  magnitudes  of  about  50-100  K  around  midnight 
[e.g.,  Faivre  et  al.,  2006],  Recent  observations  of  a  night- 
glow  brightness  wave,  commonly  associated  with  the 
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MTM,  suggest  that  the  phenomenon  may  also  extend  into 
mid-latitudes  [Colerico  et  al.,  2006], 

[3]  Mayr  et  al.  [1979]  found  a  substantial  terdiurnal 
component  in  the  AE-E  measurements,  particularly  on  days 
of  a  prominent  MTM.  They  advanced  a  mechanism  where¬ 
by  the  terdiumal  and  possibly  higher  harmonics  are  gener¬ 
ated  in  situ  by  nonlinear  interactions  between  the 
semidiurnal  tide  propagating  from  the  lower  atmosphere 
and  the  diumally  varying  ion  drag.  It  therefore  remains 
somewhat  of  a  mystery  that  subsequent  attempts  to  repro¬ 
duce  a  realistic  MTM  in  coupled  ionosphere-thermosphere 
models,  which  include  all  the  relevant  physical  processes, 
have  been  unsuccessful  even  after  substantial  adjustments 
of  amplitudes  and  phases  of  the  semidiurnal  tidal  modes 
[Fuller-Rowell,  1981;  Fesen,  1996;  Meriwether  et  al., 
2008],  It  has  been  suggested  that  higher-order  tidal  har¬ 
monics  entering  the  thermosphere  from  below  may  be 
needed  to  produce  a  realistic  MTM  in  these  models  [e.g., 
Colerico  et  al.,  2006]. 

2.  Model 

[4]  The  Whole  Atmosphere  Model  (WAM)  is  being 
developed  under  the  Integrated  Dynamics  in  the  Earth’s 
Atmosphere  (IDEA)  project  to  study  dynamical  links  be¬ 
tween  the  lower  atmosphere  and  the  upper  atmosphere  and 
ionosphere  [Akmaev  et  al.,  2008;  Akmaev  and  Juang,  2008; 
Fuller-Rowell  et  al.,  2008].  It  represents  a  neutral- 
atmosphere  component  of  the  coupled  IDEA  model  and  is 
an  extension  of  the  medium-range  weather  prediction 
Global  Forecast  System  (GFS)  model  to  the  top  of  the 
atmosphere  [Akmaev  et  al.,  2008]. 

[5]  Several  long  free-running  (i.e.,  without  constraints 
imposed  by  data  assimilation)  simulations  have  been  con¬ 
ducted  to  evaluate  the  model’s  climatology  and  internal 
variability.  The  model  was  run  for  1 5  months  and  the  output 
was  analyzed  for  the  last  12  months.  Here  we  present  a 
preliminary  analysis  of  the  MTM  found  in  the  simulations 
for  low  solar  activity  and  quiet  geomagnetic  conditions. 

3.  Results  and  Discussion 

[6]  Figure  1  presents  snapshots  of  temperature  deviation 
AT  from  the  zonal  mean  Tq  at  a  pressure  level  near  285  km 
in  March  at  20°N  (top)  and  December  at  30°S  (bottom). 
Each  thin  line  represents  one  day  in  a  given  month,  and  the 
two  columns  correspond  to  two  different  UT-longitude 
sectors.  A  prominent  MTM  with  a  magnitude  between 
about  50  K  and  well  over  100  K  clearly  stands  out  at  both 
locations  and  seasons.  Strong  day-to-day  variability  is  also 
evident,  not  unlike  the  variability  observed  from  AE-E 
[Spencer  et  al.,  1979].  An  MTM  is  invariably  seen  in  other 
seasons  at  low-  to  mid-latitudes  in  the  upper  thermosphere 
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Figure  1.  Temperature  deviation  from  the  zonal  mean  near  285  km  in  (top)  March  at  20°N  and  (bottom)  December  at 
30°S  at  (left)  UT  =  0:00  and  (right)  UT  =  6:00.  See  text  for  details. 


(for  a  sample  horizontal  structure  see  Figure  5),  sometimes 
exhibiting  more  variability  in  the  time  of  occurrence  and 
smaller  magnitudes.  In  agreement  with  the  analysis  of 
Herrero  and  Spencer  [1982],  at  a  given  height  the  magni¬ 
tude  is  typically  largest  in  summer  (e.g.,  Figure  1,  bottom), 
and  it  continues  to  increase  with  height  in  most  locations. 

[7]  A  substantial  part  of  the  internally  generated  variabil¬ 
ity  seen  in  Figure  1  is  clearly  driven  by  waves  at  relatively 
small  scales.  The  thermosphere  is  a  naturally  strongly 
dissipative  medium,  eliminating  the  need  for  “sponge 
layers”  and  extra  numerical  dissipation  often  imposed  in 
upper  layers  to  stabilize  atmospheric  model  codes.  WAM 
employs  no  “sponge  layers”  and  remains  stable  using  a 
substantially  reduced  numerical  Rayleigh  friction  coeffi¬ 
cient  even  compared  to  the  standard  version  of  GFS 
[Akmaev  et  al.,  2008],  Numerical  experiments  show  that 
both  the  size  of  the  MTM  and  especially  its  variability  are 
very  sensitive  to  imposed  numerical  dissipation.  Excessive 
numerical  dissipation  [ Sutton ,  2008]  and  perhaps  insuffi¬ 
cient  horizontal  and  vertical  resolution  may  be  some  of  the 
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Figure  2.  Temperature  amplitude  spectrum  as  a  function 
of  frequency  and  zonal  wavenumber  for  December  at  30°S 
near  285  km.  Westward  migrating  waves  are  located  on  the 
diagonal  of  the  left  plot. 
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Figure  3.  (left)  Relative  temperature  deviation  as  a  function  of  height  and  longitude  (local  time)  for  December  at  30°S  and 
UT  =  12:00.  (right)  The  same  signal  reconstructed  using  migrating  waves  with  zonal  wavenumbers  3-12  at  select  levels. 


reasons  for  the  inability  of  previous  models  to  generate  a 
realistic  MTM. 

[s]  The  whole  feature  migrates  to  the  west  following  the 
apparent  motion  of  the  sun  but  substantial  UT-longitude 
variability  is  also  evident  suggesting  a  contribution  from 
non-migrating  waves  at  various  scales.  A  double  Fourier 
analysis  in  time  and  longitude  [Akmaev  et  al. ,  2008]  reveals 
a  spectrum  consisting  of  non-migrating  waves  propagating 
both  eastward  and  westward  in  addition  to  the  migrating 
waves  (Figure  2).  Following  the  spectrum  from  the  lower 
thermosphere  shows  an  expansion  with  altitude  of  the  tail  of 
westward  migrating  waves  to  higher  frequencies  and  zonal 
wavenumbers  along  the  diagonal  in  Figure  2  (left).  These 
waves  form  the  bulk  of  the  MTM  feature  as  illustrated  by 


the  green  lines  in  Figure  1  presenting  a  monthly  mean 
reconstruction  of  the  signal  by  migrating  waves  up  to  zonal 
wavenumber  12.  The  blue  lines  represent  a  sum  of  the  first 
three  migrating  tidal  harmonics,  suggesting  that  the  terdiur- 
nal  tide  alone  is  not  sufficient  for  a  realistically  sized  MTM. 
The  red  lines  sum  up  all  waves  in  the  double  Fourier 
spectrum  up  to  zonal  wavenumber  12  and  frequency 
12  day  1  (period  2  hours).  These  are  a  good  representation 
of  the  monthly  mean  signal  and  also  reproduce  the  UT- 
longitude  variability  due  to  non-migrating  waves. 

[9]  Figure  3  presents  a  monthly  mean  relative  tempera¬ 
ture  deviation  AT/T0  as  a  function  of  altitude  and  longitude 
for  December  at  30°S  and  UT  =  12:00  (Figure  3,  left).  It 
shows  how  a  three-peak  (terdiumal)  structure  emerges  from 
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Figure  4.  Monthly  mean  temperature  amplitudes  of  migrating  waves  with  zonal  wavenumbers  2-5  as  a  function  of  height 
for  (left)  March  at  20°N  and  (right)  December  at  30°S.  The  thin  straight  lines  approximately  show  the  altitude  dependence 
for  conservative  waves. 
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Figure  5.  Monthly  mean  temperature  deviation  from  the 
zonal  mean  in  December  near  285  km  at  UT  =  12:00  as  a 
function  of  latitude  and  longitude  (local  time). 

a  predominantly  semidiurnal  variation  above  about  1 00  km. 
As  it  propagates  upward,  the  phase  (time  of  maximum)  of 
one  of  the  three  temperature  maxima  advances  from  morn¬ 
ing  hours  to  near  midnight.  This  is  also  illustrated  in 
Figure  3  (right)  presenting  a  reconstruction  of  the  signal 
by  migrating  waves  with  zonal  wavenumbers  3-12  at  select 
altitudes.  As  the  three-peak  structure  propagates  upward  and 
its  phase  advances  to  earlier  times,  two  of  the  three  maxima 
appear  to  be  suppressed  by  enhanced  dissipation  during 
daytime  and  only  the  nighttime  maximum  survives  at  a 
substantial  magnitude,  which  is  conceptually  similar  to  the 
mechanism  proposed  by  Mayr  et  al.  [1979], 

[10]  The  rapid  emergence  of  a  strong  terdiumal  tide  in  the 
lower  thermosphere  deserves  further  investigation. 
However,  it  is  already  clear  from  a  preliminary  analysis 
(Figure  4)  that  only  a  relatively  small  fraction  of  this  tide 
arrives  from  below  while  the  bulk  of  the  oscillation  is  likely 
generated  in  situ  in  the  lower  thermosphere.  This  is 
evidenced  by  the  fast  amplitude  growth  of  the  terdiumal 
wave  with  altitude  as  well  as  its  phase  discontinuities  (not 
shown)  at  the  base  of  the  thermosphere.  As  is  well  known, 
the  gravest  mode  of  the  terdiumal  tide  is  evanescent  in  the 
mesosphere  because  of  the  negative  temperature  lapse  rate 
[e.g.,  Akmaev,  2001].  In  the  thermosphere  this  mode  can 
again  propagate  freely  and  grow  exponentially  with  height 
(Figure  4).  Several  mechanisms  may  contribute  to  its 
generation  here  including  direct  solar  radiative  forcing  and 
nonlinear  interactions  between  the  diurnal  and  semidiurnal 
tides  [Akmaev,  2001],  as  well  as  possibly  nonlinear  inter¬ 
actions  between  the  semidiurnal  tide  and  ion  drag  [Mayr  et 
al.,  1979].  Relative  contributions  of  the  different  generation 
mechanisms  can  only  be  determined  by  a  further  in-depth 
analysis  of  the  respective  forcing  terms. 

[n]  Higher-order  migrating  harmonics  also  begin  to 
grow  rapidly  after  the  terdiumal  wave  attains  a  substantial 


amplitude  (Figure  4).  The  extension  of  the  spectral  tail  of 
migrating  waves  (Figure  2,  left)  most  likely  represents  a 
spectral  manifestation  of  the  selective  damping  of  the  initial 
predominantly  terdiumal  wave  during  daytime  as  it  prop¬ 
agates  upward  [e.g.,  Bracewell,  2000],  As  already  men¬ 
tioned  above,  this  is  conceptually  consistent  with  the 
mechanism  of  Mayr  et  al.  [1979]  but  applied  to  migrating 
harmonics  of  higher  order,  not  just  the  semidiurnal  tide. 

[12]  Figure  5  presents  a  horizontal  structure  of  AT  near 
285  km  in  December,  which  is  consistent  with  the  AE-E 
climatology  at  solstice  conditions  [Herrero  and  Spencer, 
1982],  including  the  common  observation  that  in  summer 
the  tropical  MTM  peaks  earlier  and  attains  larger  mag¬ 
nitudes  than  in  winter  [e.g.,  Faivre  et  al.,  2006],  As  the 
V-shape  structure  depicted  in  Figure  5  migrates  westward, 
the  MTM  will  appear  moving  poleward  to  an  observer  on 
the  ground  [Herrero  and  Spencer,  1982;  Meriwether  et  al., 
2008].  Our  simulations  also  show  the  MTM  extend  well 
into  mid-latitudes  [cf.  Colerico  et  al,  2006]  reflecting  the 
latitudinal  structure  of  the  main  terdiumal  modes  [Akmaev, 
2001], 

[13]  A  bow-wave  stmcture  of  small-scale  waves  observed 
in  nighttime  thermospheric  densities  at  400  km,  similar  to 
the  one  seen  in  Figure  5,  has  recently  prompted  Forbes  et 
al.  [2008]  to  suggest  that  these  waves  and  perhaps  the  MTM 
may  be  at  least  partially  generated  in  situ  by  the  fast 
movement  of  the  evening  terminator  through  the  thermo¬ 
sphere.  Our  simulations  clearly  demonstrate  (Figure  3)  that 
the  MTM  originates  in  the  lower  thermosphere  from  sources 
possibly  connected  to  the  middle  and  lower  atmosphere.  A 
closer  examination  of  Figure  1  also  reveals  the  presence  of 
smaller  secondary  maxima  on  some  nights  between  approx¬ 
imately  20:00-21:00  LT  as  well  as  in  the  early  morning 
hours.  Similar  features  have  been  observed  from  the  ground 
(J.  W.  Meriwether,  personal  communication,  2009).  It 
would  be  interesting  to  investigate  the  origin  of  these 
variations  and  their  possible  relation  to  the  density  waves 
observed  by  Forbes  et  al.  [2008], 

4.  Conclusions 

[14]  WAM  appears  to  be  the  first  comprehensive  model 
to  internally  generate  an  MTM  of  a  realistic  magnitude  in 
the  thermosphere.  Initial  analysis  suggests  that  the  MTM  is 
generated  in  a  chain  of  dynamical  processes  that  can  be 
traced  to  the  lower  thermosphere  and  possibly  even  further 
down  to  tidal  sources  in  the  lower  atmosphere.  WAM 
simulations  reproduce  the  salient  features  of  the  MTM 
and  exhibit  its  strong  variability.  Note  that  in  the  present 
simulations  the  variability  is  generated  internally  by  sto¬ 
chastic  tropospheric  weather  processes  alone  and  can  only 
be  expected  to  increase  with  external  forcings  modulated  by 
geomagnetic  and  solar  activity. 
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NASA  Heliophysics  Theory  and  Living  With  a  Star  (LWS)  programs  as 
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